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The diffractive production of open charm in deep-inelastic scattering is studied 
in the semiclassical approach which has been proposed recently. In this approach, 
the leading order process contains a charm quark pair and an additional gluon 
in the diffractive final state. The p^-spectrum and the diffractive mass distribu- 
tion are evaluated and compared with predictions based on perturbative two-gluon 
exchange calculations for charm quark pair production. It is shown that the p±- 
spectrum provides a clear test of the underlying partonic process whereas the 
diffractive mass distribution reflects the non-perturbative mechanism of colour 
neutralization. 



The precise measurements of the diffractive structure function at small x in the 
experiments at HERA [I] provide a challenge to the theoretical description of deep- 
inelastic scattering based on QCD. The theoretical interest, as well as the difficulty, lies 
in the interplay between 'soft' and 'hard' processes which is the characteristic feature of 
diffractive deep-inelastic scattering. To disentangle both aspects of 'hard diffraction' the 
study of less inclusive processes will be decisive. First results on the diffractive production 
of open charm have already been reported by the HI collaboration 0. 

Diffractive charm production is theoretically interesting since one might expect the 
additional hard scale, provided by the charm mass, to ensure the applicability of pertur- 
bation theory. Indeed, several authors have considered perturbative two-gluon exchange 
as a mechanism for the excitation of open charm in diffraction |3|. [|] . 

The present paper investigates charm production in the framework of the semiclas- 
sical approach to diffraction, which treats the proton as a classical colour field || ||. 
Working in the proton rest frame, the high energy scattering of a partonic fluctuation 
of the virtual photon is calculated. The leading process is the production of a ccg-final 
state, with the relatively soft gluon being responsible for the large hadronic cross section 
of the process. In our analysis we concentrate on the comparison with the perturba- 
tive two-gluon model, with a cc-final state, where the colour singlet exchange is kept 
hard by the large charm massQ This quantitative comparison assumes that higher order 
corrections to the two-gluon exchange are small for the considered observables. 

As we shall see, the p^-distribution of the produced charm jets is a clear probe of 
the underlying partonic process. The mass spectrum, on the other hand, reflects the 
additional soft partons necessarily present in the final state. These observables of diffrac- 
tive open charm production can be used to understand the importance of soft and hard 
processes in inclusive diffraction. 

Cross sections for diffractive charm production 

Recall first the qualitative picture of diffraction in the semiclassical approach |5|, || : 
for massless quarks the leading order process is the production of a colour neutral qq- 
pair off the proton colour field. In this process the leading twist contribution comes 
entirely from the aligned jet region. This can be understood intuitively, since a highly 
virtual photon can split into two nearly on-shell quarks only if they share the photon's 
longitudinal momentum in a highly asymmetric way. 

In the case of heavy quarks, the aligned jet configuration is suppressed by the quark 
mass. The reason is the large off-shellness that is necessarily present if a virtual photon 
splits into two massive particles. As a result, the production of ccg-final states, with the 
gluon being relatively soft, becomes the leading process. 

One of the two contributing diagrams is displayed in Fig. |1|. The corresponding cross 
sections, calculated in || for massless quarks, generalize straightforwardly to the massive 

1 For a general discussion and references to other approaches see @. 
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case. The results can be given in the following form, 
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Here a = Po/qo and a' = k' /q Q are the longitudinal momentum fractions carried by 
quark and gluon, Q 2 = —q 2 is the photon virtuality. All transverse momenta are defined 
relative to the 7*-direction. 




proton 



Figure |l|: Open charm production off the proton colour field. 

The a'- and /^-distributions are completely non-perturbative. They are governed by 
the Fourier transform of the colour field dependent function 

W*(yj_) = A(U\x ± + y±)U(x ± )) - 1 , (5) 

which is built from non-Abelian eikonal factors in the adjoint representation. These 
matrices correspond to light-like paths penetrating the colour field at transverse positions 
x± and x± + y±. To leading order in l/m c the two quarks have a small transverse 
separation and act like a colour octet. The two eikonal factors arise from the propagation 
of this effective octet and the gluon. 

The integration over the transverse momentum of the produced quarks is dominated 
by the hard scales Q 2 and m 2 , implying l± w —p± for most of the events. Only the 
integration region where a' and k' ± are small, in which the gluon tests large transverse 
distances in the proton, gives rise to a non-power-suppressed contribution to the total 
cross section. In this region the factor a'N 2 + k\ in the denominator of Eq. (|3|), which 
comes from the propagator of the gluon in the wavefunction of the virtual photon, be- 
comes small, ~ A 2 , where A is a typical hadronic scale. In the case where one of the charm 
quarks is soft the equivalent factor is prevented from becoming small by the quark mass, 
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one may then expand in the transverse momentum lost by the soft quark as it travels 
through the proton and the final result is power-suppressed, i.e. 0(ti 2 /m 2 ). 

The cross sections of Eqs. ([I]) and (0) can be expressed as the convolution of the 
cross section of an ordinary partonic process, namely boson-gluon fusion, with a diffrac- 
tive gluon density, which is given in terms of trW 7 ^ || . Some of the predictions discussed 
in the following analysis are not sensitive to the particular form of the diffractive par- 
ton density; in this case we refer to our approach as 'diffractive parton model'. Other 
predictions are sensitive to the eikonal approximation, which is used to treat the soft 
interaction with the proton, in this case we use 'eikonal model'. 

The following numerical analysis will compare our approach with two-gluon exchange 
calculations. The leading order cross sections for these, corresponding to the production 
of a cc-pair, are available in the literature [|3], M. Using our kinematical variables the 
results of H read 

da L 2e 2 c a em a 2 7i 2 ^G^)} 2 C n x\i n 2 ( -i (r\ 

HI - a)\ Q [a -p±) , (6) 

4(a 2 + (1 - a) 2 )p 2 ± a 4 + m 2 c (a 2 - pi) 2 } , (7) 
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where £ (or xjp) is the longitudinal momentum fraction lost by the elastically scattered 
proton. We have simplified the formulae by neglecting the scale dependence of a s and 
of the gluon density and by restricting ourselves to small transverse momenta of the 
exchanged gluons (<C a 2 +p 2 j_)- At small t the cross section is proportional to the square 
of the gluon density G(£). The factor C parameterizes the required extrapolation from 
t ~ to the integrated cross section, 
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Transverse momentum distribution and total cross sections 



The cross section formulae of the last section can be used to predict the p±- 
distributions of the produced charm jets in each model. However, in order to make a 
comparison, the integrations over the kinematical parameters of the final state gluon, a' 
and k' ± , in Eqs. (P and (@) need to be performed. 

Consider the longitudinal case, Eq. (0). Introducing the definition of W^ ± in terms 
of a Fourier transformation explicitly, the /^-integration becomes trivial. The resulting 
5-function ensures that both factors are evaluated at the same position y±. Further- 
more, it is convenient to replace the a'-integration by a u-integration, with u 2 = a'y\N 2 . 
The resulting formula is 
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where the dimensionless function 
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the M-integration is found to be logarithmically divergent. Since the eikonal approxima- 
tion will certainly fail for a gluon energy of 0(A), it is natural to introduce the cutoff 
ot' min = A/g - Assuming a soft external field the ^-integration in Eq. (||) is dominated 
by the region y\ ~ 1/A 2 , resulting in ~ x = xsj- Introducing the dimensionless 
constant 

trW x A (y A 



h A 



y± 



the leading-ln(l/x)-contribution can now be extracted from Eq. 
section follows completely analogously, and one obtains 
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These equations can be directly compared to the corresponding two-gluon results, 
Eqs. (ID and (0). In both models the a-integration can be performed analytically. 

The qualitative differences between the two models are particularly pronounced in 
the integrated cross section with a lower cut on transverse momentum of the charm 
quarks. Since most events are expected to have small y = Q 2 /(sx) (cf. 0), it is sufficient 
to consider a = o T + <j L; 

O'falmin) = [„ d P±-T2 ■ (I 5 ) 



dpi 



Fig. shows the dependence of the corresponding event fraction on the lower cut. We use 
m c = 1.5 GeV in our numerical analysis here and below. The diffractive parton model 
predicts a much stronger high-p^ tail of the distribution. For example, at Q 2 = 10 GeV 2 , 
more than 40% of the events have p\ > 5 GeV 2 , compared to only 7% in the two-gluon 
model. 

The physical reason for this difference between the two models is easily understood. 
As discussed above, the hard process in our model corresponds to boson-gluon fusion 
in the Breit frame. Therefore, the transverse momentum is logarithmically distributed 
between m 2 and Q 2 . By contrast, the two-gluon process couples the small cc-dipole 
directly to the hadron, resulting in a power suppression of small size configurations. The 
cross section comes entirely from the softest possible region, defined by the scale m 2 . No 
logarithmic tail of higher transverse momenta appears. 
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Figure 2: The fraction of diffractive charm events above p 2 ± min for Q 2 of 10 GeV 2 and 
100 GeV 2 (lower and upper curve in each pair). 

Note that our analysis does not include a s -corrections for the two-gluon exchange 
result. Such corrections have been estimated in J| and found to be important for diffrac- 
tive masses much larger than Q 2 , due to final state gluon radiation. This can affect the 
above p^-distribution of the two-gluon exchange calculation. 

For both models the complete a- and p^-integrated cross sections can be calculated 
in terms of elementary functions. Fig. |3| shows the Q 2 -dependence of the resulting ratio 
Rd = olIot- 

For our mainly qualitative discussion it is sufficient to display the Q 2 - and m 2 - 
dependence of the total cross sections in the limit m 2 <C Q 2 . In the diffractive parton 
model one finds 

°l ~ ^ , o T ~ ± (ln(Q/m c ) - |) , (16) 

while the two-gluon model results read 

<^~f(MQM)-i), °*~c?k- < 17 » 

The main qualitative difference is the suppression of the transverse cross section by 
A 2 /m 2 in the two-gluon model. Furthermore, for very large Q 2 , R£ ~ m 2 /Q 2 in the 
two-gluon model whereas in the diffractive parton model K§ ~ 1/ ln(Q/m c ). 

Charm anti-charm mass spectrum 

The different particle content of the diffractive final states, predicted by the leading 
order contributions of the two models, should be clearly visible in the resulting mass 



5 



0.25 



0.2 - 



0.15 - 



R D 



o.i - 



0.05 - ll 












two-gluon model 
diffractive parton model 



10 



40 



50 



20 30 

Q 2 (GeV 2 ) 

Figure 3: The ratio R£ = cr L /<j T for diffractive charm production. 



spectrum. In the two-gluon model the reconstructed diffractive mass should be peaked 
around the mass of the charm anti-charm pair which is given by 



M 2 = M 2 e =^±^. (18) 
a{± — a) 

In contrast, for the leading order graphs in the eikonal model, which also contain a gluon 
in the final state, the corresponding diffractive mass is given by 



k 



a 



(19) 



To quantify the above distinction between the two models, the M cg -distribution shall 
be calculated in the eikonal model for fixed M 2 . By changing variables from a' and p\ 
to M 2 and M 2 5 in Eq. (0) it is possible to derive the differential cross section 
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where a min = (1 - v / T^)/2 with 5 = Am 2 jM 2 c& and b = (M 2 - M 2 5 )/(Q 2 + M 2 e ) = 
k 2 /a'N 2 . The dimensionless function s(b) contains all the non-perturbative information 
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Figure 4: Normalized mass spectra for the eikonal model and the two-gluon model. 



and is closely related to / of Eq. 
s(b) = f d 2 k' ± (k'l) 2 



x± 



d 2 k A 
(27T) 



' 5 ii , 2k*jj b \ trW*(k' ± -k A 



k' 
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To arrive at a quantitative prediction the function s(b) has to be analyzed in more 
detail. From Eq. ([5]) and the assumption of a smooth proton colour field it can be 
concluded that trW 7 ^ 4 (y±) is a smooth localized function of y±, which vanishes together 
with its first derivative at y± = 0. This results in the limiting behaviour s(b) ~ b as 
6^0 and s(b) ~ const, as b — > oo. Since no large ratios are involved in the definition of 
s, it is natural to try the ansatz 



8(b) 



(C s + b) 



(22) 



where C s is a constant of 0(1). Fig. ^ shows a typical spectrum in M 2 5 for the eikonal 
model and illustrates that the exact value of C s is not important. The two-gluon models 
are represented by a strip at M 2 5 = M 2 . Although the width of the strip is expected to be 
of the order of the hadronic scale we have shown a larger width in the figure to illustrate 
the experimental uncertainty of the diffractive mass measurement. The normalization is 
such that the area under each curve is unity. 

A comparison of the diffractive mass, M 2 , and the charm anti-charm pair mass, M 2 E , 
for diffractive charm events at HERA should, in principle, determine which of the mech- 
anisms is responsible for the bulk of the events. 
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Ratio of charm in diffraction 



Finally, we estimate the fraction of charmed events in the total diffractive cross 
section. Since this will involve strong assumptions about the soft part of the eikonal 
model, the results of this section are less reliable than the previous part of the paper. 

The leading order diffractive cross section, determined by the production of a pair 
of light quarks in an aligned-jet type configuration, reads M 



do T _ a em 2 
da'dk? ~^7^ 6q 



d 2 k ± k_ L tiW£ L (k , ± — kx) 



(2vr) 2 a'Q 2 + k] 



(23) 



where k' ± and a' are transverse momentum and longitudinal momentum fraction of the 
outgoing soft quark and the sum is over the light quarks u, d, s. The function W T is 
defined as in Eq. (H), but with the colour matrices in the fundamental representation. 

The corresponding total cross section, calculated as above, is found to be 

where the constant hjr is defined analogously to Eq, 

This has to be compared to the charm cross section, obtained by integrating Eq. ([14]). 
To leading ln(l/x) ln(Q 2 )-accuracy the result reads 

ccg & c CX-em a s ln(l/x) . . . /r»e\ 

° T = ^sg2 \n(Q/m c )h A . (25) 

The appearance of the above large logarithms suggests that an analogous qqg-cross 
section for light quarks will form an important correction of Eq. ( |24j ) . However, this cross 
section is plagued by an infrared divergence in the region where all three particles have 
small p±. The leading-log contribution calculated in can also be obtained by replacing 
the heavy quark mass in Eq. ( p5f) with an infrared cutoff equal to Aqcd, 

afg = a em a s Jn(l/x) ^ ^ ^g/^^ . (26) 

Neglecting the small longitudinal contribution (cf. Fig. |^) the ratio r of charm in 
diffraction can now be given, 



U rp 



a f + a f 9 + a c T C9 ' 



(27) 



A numerical evaluation within the present model requires the ratio hjs./hjr. To estimate 
this ratio it is convenient to introduce the dimensionless functions 



jn(ylA 2 ) = A 2 [ trW£(<a) 2 , (2 

J X 
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where TZ G {A, J 7 } labels the representation. For smooth colour fields of finite transverse 
extension Eq. @ implies jn(z) ~ a\z 2 for z <C 1 and jn(z) ~ 6^ for 2 ^> 1, with two 
unknown constants a-ji and b-R. 

In the region of small z the matrices {7 and W in the definition of W can be expanded 
in powers of the gauge field G^. The first term contributing to tnW is of order z 2 G 2 , and 
standard formulae of representation theory |§ give a A /ajr = 2N for the group SU(N). 
In the region of large z, corresponding to large y in Eq. (|5|), it is natural to assume that 
the matrices U and U* are not correlated. Note that in Eq. ([5]) an averaging procedure 
over all colour field configurations of the proton is understood || . In the case of strong 
fields the simplest assumption is that U and W are uniformly distributed over SU(N), in 
which case the first term of Eq. (|5|) vanishes under the above average. The second term is 
simply proportional to the dimension of the representation, so that b A /bjr = (N 2 — 1)/N. 

Although the constants h-ji are given by the simple formula 

hn = nj -^jn(z) , (29) 

knowing the limiting behaviour of j A {z) an d jr{ z ) does not yet imply a knowledge of 
hyx/hjr- The problem is that the form of the interpolating function influences, in general, 
the integral in Eq. fl29|) . Assuming, for simplicity, that the functional form does not 
depend on the representation, i.e. 

3A{z) = {b A jb T f 3T{z{a A b T )j{a T b A )) , (30) 

the integrals for TZ = A and 1Z = T are found to be related by a multiplicative factor, 
h A /hr=(a A b A )/(a r b r ) = 2{N 2 -l). 

Specifying iV = 3 and using x = 10" 3 , Q 2 = 36 GeV 2 , a s (Q 2 ) = 0.22, and A QCD = 
200 MeV, the formally leading cross section in the massless quark case is found to be 
small, a^/aj? 9 m 0.1. Therefore the diffractive production of massless quarks is dom- 
inated by the qqg component of the 7* wavefunction. It is then natural to expect the 
ratio of charm in diffraction to be similar to inclusive deep-inelastic scattering (cf. |T0|). 



Indeed, explicit evaluation of Eq. (\Z7\), with the parameters given above, results in a 
charm fraction of r w 0.2. 

Conclusions 



The p_i_-spectrum and the diffractive mass distribution for the production of open 
charm have been evaluated in the semiclassical approach. The resulting cross sections 
can be expressed as the convolution of a cross section of an ordinary partonic process, 
namely boson-gluon fusion, with a diffractive parton density. In this way the physical 
content of the results becomes most transparent. 

The p^-spectrum provides a clear test of the underlying hard partonic process. It 
is found to have sizeable contributions from the region p 2 ± ~ Q 2 , which are absent at 
leading order in the two-gluon model. 
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The diffractive mass distribution, on the other hand, reflects the soft interaction of 
the proton, which is treated in the eikonal approximation. The presence of the final state 
gluon yields a production cross section not suppressed by A 2 /m 2 and a large ratio of 
charm in diffraction, comparable with the corresponding ratio in inclusive deep-inelastic 
scattering. It also allows a large difference between the diffractive mass M 2 and the invari- 
ant mass of the charm jets M 2 g . In contrast, both masses are equal, up to hadronization 
effects, for pure cc-final states. 

In summary, several observables have been identified which can discriminate between 
soft and hard mechanisms for colour singlet exchange in diffractive charm production. 
Diffraction with high-p^ jets is kinematically very similar, since the hard scale m c is 
merely replaced by p±. Therefore, we expect that the study of diffractive charm pro- 
duction will help to clarify the relative importance of soft and hard contributions to 
diffraction in general. 

We would like to thank J. Bartels, A.D. Martin, H.P. Shanahan and T. Teubner for 
valuable discussions and useful comments. 
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